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motor detachment and attachment, in steady state of the system, this three-state 
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phases) and boundary layers. In this study, using mean-field analysis, the conditions 
of existence of shock and boundary layers are obtained theoretically. Combined with 
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helpful to understand the biophysical properties of the collective transport of kinesin 
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I. INTRODUCTION 



Molecular motors are biogenic force generators acting in the nanometer range. 
They are responsible for intracellular transport of wide varieties of cargo from one 



location to another in eukaryotic cells 



-|6|. Linear motors produce sliding move- 



ments along filamentous structures called protein tracks; for 



along actin filament 
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10] , kinesin f lMl5| and dynein 16M20] along microtubule. 



example, myosin slides 



Microtubule and filamentary protein actin are protein filaments which form part of 
a dual-purpose scaffolding called cytoskeleton, act like struts or girders for the cel- 
lular architecture and, at the same time, also serve as tracks for the intracellular 
transportation networks. 

However, experiments found, one single filamentary track is usually traveled along 
by multiple motors. Fundamental understanding of these collective physical phenom- 
ena may expose the causes of motor-related diseases (e.g., Alzheimer's disease). In 



literatures, these phenomena are usually described by the totai 



exclusion process (TASEP), which is originally proposed in 21], consists of particles 



ly asymmetric simple 



hopping unidirectionally with hard-core exclusion along a ID lattice. TASEP is one 
of many examples for driven systems with stationary nonequilibrium states, which 



cannot be described in terms of Boltzmann weights. To mode 



the attachment and 



detachment of motors to and from tracks, Parmeggiani et al [221 . 123] discussed a class 
of driven lattice gas obtained by coupling ID TASEP to Langmuir kinetics, in which 
the attachment and detachment of motors is mode 
hilation respectively. Furthermore, Lipowsky et al 
model, in which the diffusion of motors in the cell is considered explicitly. However, 
in reality, a motor protein is not a mere particle, but an enzyme whose mechanical 
movement is coupled with its biochemical cycle. Therefore, recently, a three-state 



ed as particle creation and anni- 
24 . |25] suggested a more general 



model is presented by Nisi 
headed kinesin KIF1A 1 26 



linari et al to describe the intracellular transport of single- 
'I 

27] . In which, the microtubule (MT) binding motor might 
be in two states: strongly MT binding state and weakly MT binding state, denoted 
by S, W respectively. Biochemically, the strongly binding state corresponds to bare 
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motor or ATP binding motor state, and the weakly binding state corresponds to ADP 
binding state. 

One of the important feature of the collective motion of motors along one single 
track is the possible appearance of shock or domain wall, which is defined as the 
interface between the low-density and high density regions. For the usual two-state 
TASEP, using mean field method, the existence and properties of the shock have 



been discussed recently 



28] . In this study, similar analysis to the Nishinari's three- 



state model will be presented, the efficient and necessary conditions of the existence 
of shock will be given theoretically, and with the aid of numerical calculations the 
properties of the shock will also be discussed. The method used in this study can be 



28|. 



regarded as a generalization of the one presented in 

In the next section, the three-state model and its mean field approximation will 
be briefly introduced, and then in Sec. Ill the conditions of the existence of shock 
will be presented. The existence of boundary layers and the properties of shock will 
be discussed in Sec. IV and V. Finally, this study will be shortly summarized in Sec. 
VI. 



II. TASEP WITH THREE INTERNAL STATES 

The three-state TASEP given in 



27] can be mathematically desccribed as follows. 
Let Si and Wi denote the probabilities of finding a molecular motor in the states 1 
and 2 at the lattice site i at time t respectively (states 1 and 2 correspond to strongly 
bound and weakly bound states of molecular motors). Then Si, Wi are governed by 
the following master equations 

— =u a {l - Si - Wi) - u h Si - u d Si + u s Wi 
dt 

+ UfWi-til -Si- Wi) + (1 - c)u f Wi{S i+1 + Wi+i), 
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^ = - (u a + w f )Wi(l - S i+1 - W l+1 ) + uiJSi 
-[u s + (l-c)u f ]Wi(S i+1 + W i+1 ) 

(2) 

- u b Wi(2 - S i+1 - W i+1 - SU - Wi-x) 

+ u b (W i - 1 + W i+1 )(l - Si-Wi), 
where u a is the rate of a molecular motor binding to the empty lattice site i, i.e. the 
transition rate of state to state 1, Uh is the transition rate of state 1 to state 2, i.e. 
the rate of ATP hydrolysis, Ud is the transition rate of state 1 to state 0, i.e. the 
rate of detachment, u b is the rate of random Brownian motion. After the release of 
ADP, the motor steps forward to the next binding site in front with rate Uf, stays at 
the current location with rate u s . c is an interpolating parameter (0 < c < 1). The 
corresponding equations for the left boundary (i = 1) are given by 

' ' 1 a (l-3 1 -W 1 )-u h 3 1 -'y 1 3 1 +u s Wi + {l-c)u;fW 1 (3 2 + W2), (3) 



dt 

dW x 
~df 



+ u f )W 1 + + ow / Wi(5 , 2 + W 2 ) 



(4) 



+ u b W 2 (l -St- Wi) - u b Wi(l -S 2 - W 2 ) - l2 W x , 
where a is the rate of attachment of the motors at the left boundary (i.e. the lattice 
site i — 1), 7i, 72 are the rates of detachment of motors in state 1 and state 2 at the 
left boundary respectively. The equations for the right boundary (i = N) are give by 
dSpj 



(5) 



(6) 



(H --5(1 -S N - W N ) + u f W N - X {l -S N - W N ) 
+ lu s W n - u h S N - PiSn, 

dW N 

—— =u h S N - uj s W n - f3 2 W N + w 6 Wat_i(1 - S N - W L ) 

- u b W N {l - SV-i - W N -i), 
where 5 is the rate of attachment of the motors at the right boundary (i.e. the lattice 
site i = N), Pi (3 2 are the rates of detachment of motors in state 1 and state 2 at the 
right boundary respectively. 

It should be pointed out that the exclusion process described above is different from 
the one discussed in [29] , where multiple occupancy of sites is allowed if particles are 
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in different internal states. Here, the multiple occupancy is unallowed. However, 
the particles bounding to the lattice site might be in two different states 1 and 2, 
corresponding to the strongly bound and weakly bound states. As in 22] , attachment 
and detachment of a motor are modeled as, effectively, creation and annihilation of 
motors on the lattice. Moreover, the transition between states 1 and 2 is described 
by the rates Uh, oo s , Uf, the Brownian ratched mechanism is described by rate Ub- 

Mean Field Approximation: In the large N limit, we can make the continuum 
mean field approximation to Eqs. and (J2]). Let Ax = -^-j- and x = (i — I) Ax. 
Obviously, < x < 1, since 1 < i < N . Using the Taylor expansion 

k=l ' k=l 



(7) 



The continuum limits of Eqs. ([I]) and ([2]) are then 

=w„(l - S - W) + a.W - K + w)s 

+ (1 - c).,W (S + W + Ax^l + Ax^M 
+ 0(Ax 2 ), 

dW(x,t) „ / . <9S(x,t) . dW(z,t)\ 

=cw/W (5 + W + Ax — y-i- + Ax y^ 

at \ ox ox J 



(9) 



- (u s + tj f )W + u h S + 0(Ax 2 ). 
Thus, the probability density p(x, t) = S(x, t) + W(x, t) of finding a molecular motor 
at lattice site x at time t satisfies [summing Eqs. (jSD and 0] 

^^(l-,)-,^ + u ,»5^A« + 0(A*»). (10) 



As the discussion in [30|, in the thermodynamic limit N — > oo, there are three 
regimes to be distinguished. If u a and are of order [1/(N — l)] a with a < 1, then 
at the steady state, the system, Eqs. (jHJ) and f lTOj) . reduces to 

^(1-^-5) -^ = 0, 

(11) 

cw/H/ (S + W)- (cu s + u f )W + co h S = 0. 
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So the probability S satisfies 

ck{k + 1)S 2 + [(Jfe + + 1) + n - c(2k + 1)}S + (c - I - 1) = 0, (12) 

and W — 1 — (k + 1)S, p — S + W, where I = u s /uf, n = Uh/uf, k = Ud/oj a (because 
of the particle-hole symmetry, we restrict the discussion to the case u a > Ud, i.e. 
< k < 1). 

For the cases that u a and Ud are of order [1/(N — 1)] Q but with a > 1, the local 
kinetics is negligible and the system will be 

W(W + S-1) =c, 

(13) 

cuj f W (S + W) - (co s + u f )W + u h S = 



where the constant C is determined by the left or right boundary conditions [3l|, |32|. 
The case of the local rates u a and Ud being of the order 1/(N — 1) is the most 
interesting one, and will be investigated further in this study. In the following, we 
always assume that the local rates u a and Ud are of the order 1/(N — 1). 



Let 



III. THE EXISTENCE OF SHOCK 



= ^= ^, (14) 

Ax' Ax' v ; 



then at steady state, the leading terms of Ax of Eqs. OH]) flHIj) are 

coo f W (S + W)- (u s + u f )W + u h S = 
or 

u f dw{p d ~ l) + n a (i - P ) - n d ( P -w) = o, 

cu f Wp - (uj s + uj f )W + u h (p -W) = Q. 
The second equation implies 



(15) 



(16) 
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The steady state flux is then proportional to 



J = W{1- p) 



np{l - p) 



n + I + 1 — cp 

n 

The same as in [271 ]. at the steady state, we can obtain the left boundary conditions 

a — [ca(a — u s )/uf] 



5(0) 

and the right boundary conditions 



W(0) 



S(l) 



0J s + /3 



W(l) 



a: 

UJ f '' 



(19) 



/3 



(20) 



u h + u s + P UJf' 

One can find that Eqs. ( JT6l) involves only the first-order derivatives of p and 
with respect to x whereas there are two sets of boundary conditions (|T9|) and 
( 1201) . Therefore, if we integrate the equations ( JT6l) with the left boundary conditions 
( TT9l . the solution (denoted by pi, Wi, Si respectively) may not, in general, match 
smoothly with the solution (denoted by p r , W r , S r respectively) obtained for the 
same equations but with the right boundary conditions ( !20l) . This discontinuity 
corresponds to a shock or domain wall. However, at any position x, the continuity 
condition of motor flux, or equivalently J/(x) = J r (x), should be satisfied, where 
Ji(x) = W(x— )[1 — p(x— )] and J r (x) = W(x+)[1 — p(x+)]. At the shock position x s , 



W(x s -) = Wi(x s ), p{x s -) = pi(x s ), 

W(X S + ) = W r (x s ), p{x s + ) = p r {x s ). 

So the continuity condition Ji(x s ) = J r (x s ) implies 

npi(x s )(l - pi{x s )) N _ np r (x s )(l - p r (x s )) 

■ J [Xg 



(21) 



n + l + l — cpi(x s 



n + l + l — cp r (x s 



or 



where 7 



n+l+l 



Pl(Xs) + Pr(Xs) = 1 + -J(X S ) = l + jpi(x s )p r (x s ), 

n 

< 1. 



(22) 



(23) 



From Eqs. (fl6l) . one can show the probability p satisfies 



(7P 2 - 2p + l)p x = n ah (l - 7p) 



(k + l)cp 2 - - + c + k(l + 1) )p + - 

7 / 7. 



(24) 



s 



with Q a h = Q a /u>h- It can be proved that, for < c < 1 and I > 0, the discriminant 



A := 



- + c + k(l + l) 

7 



4{k + l)c 2 

7 



> 0. 



So the equation ( 1241) can be reformulated as 



(p - Ps)(p - P4)Px = -fiah(A; + l)c(p - p )(p - Pi)(p - P2), 



(25) 



where 



Pl,2 



7 



P3,4 



1 

7 



^ + c + fc(/ + l) T W £ + c + A;(/ + l) 



4(fc+l)c 2 

7 



2(Jfe + l)c 

One can easily show that, for < c < 1, / > and k > 0, 



1 

P0 > 1, P4 > 1, P2 > 1, Pi < 1 ^ - ^ - L 



Particularly, 



Pi 



2c 

7 



+ c + /c(Z + l) 



1 



Z+n+1 



+ c + /c(Z + l) 



as k — t- 0, 



4(fc+l)c 2 



< 1 as c -»■ 0. 



(26) 



(27) 



(2c 



Z+n+l+fc(i+l) 

Moreover, one can easily show that the function f(x) = x — y/ x 2 — x decreases with 
x > 1 monotonously. So, for c < (1 — k) (I + 1) + n, 



Pi 



± + c + k(l + l) 



± + c + k(l + l) 



4(fc+l)c 2 

7 



2(& + l)c 



(k + l)(l + l)+n + c f(k + l)(l + l)+n + c\ 2 l + n + 1 



2(k + l)c 
{k + l)(l + l)+n + c 



2(k + l)c 



(k + l)c 



> 



{k + l){l + l)+n + c\ 2 (k + l){l + l)+n + c 



> 



2(k + l)c 

l + n + 1 

c 

=P3- 



2(k + l)c 



2(& + l)c 



Z + n + 1\ l + n + 1 



(29) 
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Therefore, in the following, we always assume that p , p 2 , Pa > 1, and ~ < p 3 < p\ < 1. 
The general solutions of Eq. ( l25l) are 

F(p) =x + C, (30) 



where C is an arbitrary constant and 

1 

n ah (k + i)c 



F (P) = - o //, r\ l A ln \P - Pol + B ln IP ~ Pil + D In \p - p 2 \] , (31) 



with 



A = Po(Po - 1) B = P\ ~ 2poPi + Po ^ _ p| - 2p P2 + Po 



(PO - Pl)(P2 - PqY (PO -Pl)(Pl - Pi)' (PO -P2)(P1 -P 2 )' 

(32) 

So the solution of Eq. f l23|) . which satisfies the left boundary condition p(0) = 5(0) + 
W(0), see Eq. ([19]), is 

F( Pi ) = a; + Fo, or Pl {x) = F~\x + F ), (33) 

where F = F[p(0)]. Similarly, the solution of Eq. (l2"o"|) . which satisfies the right 
boundary condition p(l) = S(l) + 1^(1), see Eq. (120]) . is 

F(p r ) = x + F x - 1, or p r (x) =F- 1 (x + F 1 -l), (34) 

where iq = F(p(l)). In the following, we assume p\ ^ p r (otherwise, there would be 
no shock and boundary layers). 

Combining ( 121]) ( 123]) (1331 (|34j) . one sees that, at the shock position x s 

F( Pl (x s )) - F (p r (x s )) +Fx-F o -l = 0, 

(35) 

pi(x s ) + p r (a; s ) = 1 + -ipi(x s )p r (x s ). 
These are the efficient and necessary condition of the existence of shock at the position 
x s . In other words, at the shock position x s , H(x s , fl a , k, uj, I, n, c) := 

F( Pl (x s )) - F ( !~~ Pl ^ s \ ) +Fx-F o -l = 0. (36) 



1 - -ypi(x s ). 

If there exists < x s < 1, such that H(x s , Q a , k, Uf, I, n, c) = 0, the shock will appear 
at x s , and the height of the shock is 

1 - 2pi(x s ) +jpf(x s ) 



\pr(x s ) - Pl{x s )\ 



l-jpi(x s ) 



(37) 
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Let 

H{p) := F{p) - F 

then from (13TI) . one can obtain 
dH{p) _ dF(p) _ d ( 1 - p 



P 



1-7P 
7 p 2 - 2p + 1 



F x - F - 1 = 0, 



(38) 



dp dp dp* \l->ypj n ah (l-^yp) 

[(1 - p) - Pl (l - 7 p)][(l ~ p) ~ p 2 (l ~ 7P)] + ( 7 ~ 1)(P - Pi)(p - p 2 ) 
(p ~ Pi)(P ~ P2)[(l - p) ~ Pi(l ~ 7P)][(1 - P) ~ P2(l - 7P)] 



7(l-p) 2 (l-7Pi)(p-P3)(p-p 4 ) i 



i-pi 



-7Pi 



- P2 



- 7P) 3 (P - Pi)(p " P 2 ) - Pi) (y^ - P 2 ) 



> 
< 

> 



if < p < £ 



i-pi 



if 



i-pi 



7Pi ' 

1_ 7P1 < P < Pi, 

if Pi < P < 1- 



(39) 



Using this property of the function H(p), we can obtain the following results: 

(I) For < p(0) < p 3 and p 3 < p(l) < p 1; the conditions of existence of shock in 

interval (0, 1) is 



1 



pan 



< F + 1 and F 



P(0) \ 



<F X -1, 



(40) 



i- 7 p(i)y Vl-7P(0)7 

see Fig. [T] (a). From (I25p (13~Tj) . one can find the function F(p) increases with p for 
< p < P3 and px < p < 1, and decreases with p for p$ < p < p\. Thus 



F 
F 



1-P(1) 
1-7P(1) 
1 ~ P(0) \ 



<F + 1 
< Fi - 1 



Pi 



Pr(0) < 



l- 7 p(l) 
1 - p(0) 



< 1, 



(41) 



1-7P(0); ~* ~ ' ' rrvv ^ 1-7P(0)' 
Therefore, the conditions presented in ( |4"0]) are generalizations of the ones obtained 
in 281 for the usual TASEP with motor detachment and attachment: 



pr 1 (l-p(l)) < 1, and p r (0)<l-p(0). 
In fact, if the parameter c = (i.e. 7 = 0), the Eq. ff24|) reduces to 



(42) 



1 - 2p)p x = tt ah [{l + n + l)-[{l + n + l)+k{l+ l)]p], 



(43) 
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which is similar as the model discussed in [22[ for the usual TASEP. For such reduced 
cases, P3 = 0.5, and the conditions ( 140]) of existence of shock is reduced to ( )42l) . 
(II) For < p(0) < p3 and p\ < p(l) < 1, the conditions of existence of shock in 
interval (0, 1) is 

see Fig. [1] (b). Similar as in (I), 

Therefore, conditions (I44p are also generalizations of the ones for the usual TASEP 



(45) 



pr 1 (l-p(l))<l, and Pr - 1 (l-p(0))>0. (46) 

(III) For < p(0),p(l) < p3, the condition of the existence of shock in (0, 1) is 

Pr 1 (p 3 )<l and p r (0) < \Z^$) (47) 
where p r {x) is one of the solutions of differential equation ( 124"|) . which satisfies p r (l) = 
p 3 and p r (0) > p 3 . See Fig. [1] (c). 

(IV) For p(0), p(l) > p3, there exists no shock in (0, 1). It can be readily verified that 
the function f(pi,p r ) = Pi + Pr — JPiPr — l increases monotonously with < pi,p r < 1, 
and /(p3,p3) = 0. For p(0),p(l) > P3, one knows that p 3 < min(p/,p r ) < max(p/,p r ) 
(note: we always assume pi ^ p r ). Thus f(pi,p r ) = Pi + Pr — IPiPr — 1 > 0. It is to 
say that there exists no shock [see (J33J]. 

(V) For < p(l) < p 3 and pz(0) > p 3 , there exists no shock in (0, 1). 

In conclusion, p(0) < p3 is one necessary condition of the existence of shock in 
(0,1). 



IV. THE EXISTENCE OF BOUNDARY LAYER 



Generally speaking, if p\ ^ p r and there is no shock in (0, 1), boundary layer will 
appear at least at one of the boundaries and 1. Similar as in [28[, we have the 
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following results: 

(I) For < p(0) < p 3 and p 3 < p(l) < Pi- if 



pan 



>F + 1 and f(1^)< Fi -1, (« 



1-7P(1)7 

there exists boundary layer at the right boundary x = 1, see Fig. [T] (d); if 

Vi-7P(i)y \i-ip(o)J 



(49) 



there exists boundary layer at the left boundary x = (in these cases, the shock 
position x s < 0), see Fig. [2] (b). 

In view of the property (1391) of function H(p), the conditions ( HHl) can be simplified 

as 

'>F„+1, (50) 



.1-7P(1) 

and the conditions f|451 can be simplified as 



(II) For < p(0) < p 3 and pi < p(l) < 1: if 

there exists boundary layer at the left boundary x = (i.e. the shock position x s < 0), 
see Fig. [2](c); if 

Kl^T)H + 1 ' (53) 

there exists boundary layer at the right boundary x = 1 (i.e. the shock position 
x s > 1), see Fig. [2] (a). 

(III) For < p(0),p(l) < p 3 and pz(l) ^ p(l), there exists boundary layer at x — 1. 
See Fig. GO (c), Fig. [2] (d) and Fig. |3] (a). If 

pr 1 (ps)<l and p r (0)> ^J^ , (54) 
there is also the boundary layer at the left boundary x = 0. 
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(IV) For p 3 < p(0) < 1,0 < p(l) < p 3 and p(0) 7^ p r (0), there exist boundary layers 
at both x = and x = 1, see Fig. [3] (b) and Fig. |3] (c). For these cases, p(x) = p r (x) 
for < x < 1. 

(V) For p 3 < p(0),p(l) < 1 and p(0) 7^ p r (0), there exists boundary layer at x — 0. 
For these cases, p(x) = p r (x) for < x < 1, see Fig. [3] (d) and Fig. HI 



V. THE PROPERTIES OF SHOCK 



Finally, we discuss the properties of shock briefly. From the discussion in Sec. Ill, 
one knows that p(0) < p 3 is necessary for the existence of shock. So, at the shock 
position x s , pi(x s ) < p r (x s ), see Eq. (123]) . and the height of the shock is [see Eq. ([3j 



p r {x a ) ~ pi{x s ) = ^ f V ■ ( 55 ) 



The derivative of the height e s with respect to pi(x s ) is 

de s 7-I 



dpiixg) (1 - <ypi(x s )y 



1 < 0. (56) 



At the same time, p t (0) = p(0) < p 3 means > [see (B51)]. So 



de s de s dpi{x s 



< 0, (57) 



dx s dpi(x s ) dx & 

which means, the shock height e s decreases with the shock position x s . Therefore, we 
only need to give the relations between shock position x s and the model parameters 

UJf, U s , Uh, Wft, c. 

Because of the complexity of the function F [see QQ], it is difficult to get theoteti- 
cal results as in [28]. From numerical calculations, we find that, the shock position x s 
decreases with parameters Q a , a, Ub, oj s , but increases with parameters f^, (3, cjf, Uh, c, 
see Fig. |5j In the calculations, x s is obtained by (|36|) and (|33|) . 



VI. CONCLUDING REMARKS 

In this study, the three-state process, which is presented in {26I, |27] to luodcl 
the intracellular transport of single-headed kinesin KIF1A , is theoretically analyzed 
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using mean field approximation. By similar methods as for the usual TASEP [28|, the 
conditions of the existence of shock or domain wall, which is defined as the interface 
of low-density and high-density phases, are obtained. With the aid of numerical 
calculations, the parameters dependent properties of the shock are also discussed. 
The results obtained in this study will be helpful to understand the real biophysical 
properties of motor traffic in eukaryotic cells. 
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